Lifetime of dynamical heterogeneity in a highly supercooled liquid 
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We numerically examine dynamical heterogeneity in a highly supercooled three-dimensional liquid 
via molecular-dynamics simulations. To define the local dynamics, we consider two time intervals, 
Ta and Tngp. Ta IS the a relaxation time, and Tngp is the time at which non-Gaussian parameter of 
the van Hove self-correlation function is maximized. We determine the lifetimes of the heterogeneous 
dynamics in these two different time intervals, rhetero(T"a) and rhctero(Tngp), by calculating the time 
correlation function of the particle dynamics, i.e., the four-point correlation function. We find that 
the difference between rhetero(T"a) and rhotGro(''"ngp) increases with decreasing temperature. At low 
temperatures, rhctcro(T"a) is considerably larger than Ta, while rhctcro(T"ngp) remains comparable to 
Ta- Thus, the lifetime of the heterogeneous dynamics depends strongly on the time interval. 



One of the long-unresolved problems in material sci- 
ence is the glass transition [1] [2]. In spite of the ex- 
tremely widespread use of glass in industry, the forma- 
tion process and dynamical properties of this material 
are still poorly understood. Numerous studies have at- 
tempted to explain the fundamental mechanisms of the 
slowing of the dynamics observed in fragile glass (i.e., 
the sharp increase in viscosity in the vicinity of the glass 
transition). However, the physical mechanisms behind 
the glass transition have not been successfully identified. 

Recently, dynamical heterogeneities in glass-forming 
liquids have attracted much attention. In a system dis- 
playing dynamical heterogeneity, the dynamic character- 
istics (i.e., particle displacements and local structural 
relaxations) are non-uniformly distributed throughout 
space. Dynamical heterogeneities have been detected and 
visualized through simulations of soft-sphere systems [3] 
[4] [5] [6] [7] [8] , hard-sphere systems [9] , Lennard- Jones 
(LJ) systems [10], and experiments [11] [12]. Insight into 
the mechanisms of dynamical heterogeneities will lead to 
a better understanding of the slowing of the dynamics 
near the glass transition. 

Conventional two-point density correlation functions 
are not informative when applied to the investigation 
of dynamical heterogeneities. We need to examine the 
correlation of the particle dynamics, not just snapshots. 
We can quantify the correlation length ^ of the hetero- 
geneous dynamics by calculating the four-point correla- 
tion functions, which correspond to the static structure 
factor of the particle dynamics. Several simulations [5] 
[13] [14] [15], experiments [16] [17], and mode-coupling 
theory [18] have estimated ^ in terms of the four-point 
correlation functions and revealed that ^ increases with 
decreasing temperature. In addition, we can quantify the 
lifetime Thctoro of the heterogeneous dynamics by employ- 
ing the multiple time extension of the four-point correla- 
tion functions (i.e., the multi-time correlation functions), 
which correspond to the time correlation functions of the 



particle dynamics. Thetero has been measured in terms of 
the multi-time correlation functions by simulations [5] [6] 
[19] [20] and experiments [16] [21] [22]. It was reported 
that Thctcro increases dramatically with decreasing tem- 
perature and can become greater than the a relaxation 
time near the glass transition. 

In 2009, Kim et al.[20] investigated the correlations 
between the heterogeneous dynamics at various time in- 
tervals. They calculated the sum of the time correla- 
tion functions for the heterogeneous dynamics at these 
time intervals and determined the lifetime of the hetero- 
geneous dynamics as a characteristic time at which the 
sum of correlation functions decays. 

In this letter, we demonstrate via molecular-dynamics 
(MD) simulations that the lifetime of the heterogeneous 
dynamics depends strongly on the time interval. To de- 
fine the local dynamics, we consider two time intervals: 
the a relaxation time, r^, and the time Tngp at which the 
non-Gaussian parameter of the van Hove self-correlation 
function is maximized. We estimate the lifetimes of the 
heterogeneous dynamics in these two different time inter- 
vals ,rhotoro(TQ) and Thctcro (^ngp), by Calculating the time 
correlation function of the particle dynamics. Finally, we 
compare the two lifetimes. 

The conventional two-point correlation function 
F{k, t) represents the correlation of the local fluctuations 
Sn{k,t) in some order parameter, such as the particle 
density. 5n{k,t) is the Fourier component k of the fluc- 
tuations at the time t, and F{k,t) = {dn{k,t)6n{—k,0)), 
k — \k\. The two-point correlation function can describe 
the particle dynamics in the time interval [0,t\, averaged 
over the initial time and space. As the time interval t in- 
creases, F{k, t) decays in the stretched exponential form. 
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where T(fc) is the relaxation time of the two-point cor- 
relation function, which represents the characteristic 
timescale of the averaged particle dynamics. To exam- 
ine the lifetime of spatially heterogeneous dynamics, we 
have to calculate the time correlation function of the 
local fluctuations SQk{q,to,t) in the particle dynamics. 
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FIG. 1. Schematic illustration of two time intervals and their 
time separation. 



SQk{q,tQ,t) is the Fourier component q of the fluctu- 
ations in the particle dynamics associated with a mi- 
croscopic wavenumber k in the time interval [to, to + t]. 
F{k,t) is equal to Qk{q,to,t) averaged over the initial 
time to and space, i.e., F{k,t) ^ {Qk{q,to,t)). The time 
correlation function defined by 

Fi^kiq, ts,t)^ {SQk{q, ts + t, t)dQk{-q, 0, t)), (2) 

represents the correlation of particle dynamics between 
the two time intervals [0, t] and [tg +t,ts+ 2t]. ts is the 
time separation between the two time intervals [0, t] and 
[ts + t,ts + 2t], as is schematically illustrated in Fig.l. 
Fi,k{<l,ts-,t) is the multiple time extension of the four- 
point correlation function [20] . As the time separation ts 
increases, i^4^fc(g, is, t) with fixed t decays in the stretched 
exponential form. 
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where T4_kiq,t) is the relaxation time of the correlation 
of the particle dynamics. We determined the lifetime 
of the heterogeneous dynamics Thetero(0 '''4,kiq,t) at 
q — 0.38, the smallest wavenumber in our simulation. 

To calculate the time correlation function of the par- 
ticle dynamics, we performed MD simulations in three 
dimensions on binary mixtures of two different atomic 
species, 1 and 2, with Ni = N2 = 5,000 particles and a 
cube of constant volume V as the basic cell, surrounded 
by periodic boundary image cells. The particles inter- 
act via the soft-sphere potentials Vabir) = ^{(JablrY'^, 
where r is the distance between two particles, tJah = 
((Ta + crft)/2, and a,b d 1,2. The interaction was trun- 
cated at r = ScTafc. In the present letter, the following 
dimensionless units were used: length, cti; temperature, 
e/ks] and time, tq = {mia\/e). The mass ratio was 
m2/mi = 2, and the diameter ratio was a^jai — 1.2. 
This diameter ratio avoided system crystallization and 
ensured that an amorphous supercooled state formed at 
low temperatures [23]. The particle density was fixed at 
the high value of p = (iVi + N2)(y\ly = 0.8. The system 
length was L = = 23.2(71. Simulations were carried 
out at T = 0.772, 0.473, 0.352, 0.306, 0.289, 0.267, and 
0.253. Note that the freezing point of the corresponding 
one-component model is around T = 0.772 (Feff = 1.15) 
[23]. Here, Tqs is the effective density, a single parameter 
characterizing this model. At T = 0.253 (Fes = 1.52), 
the system is in a highly supercooled state. We used the 
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FIG. 2. Ta, Tngp versus inverse temperature l/T. We use 
these two time intervals to define the local dynamics. 



leapfrog algorithm with time steps of 0.005t when inte- 
grating the Newtonian equation of motion. Very long 
annealing times (3 x 10^ for T = 0.253) were chosen. No 
appreciable aging effect was detected in various quan- 
tities, including the pressure or the density correlation 
function. 

To define the local dynamics, we consider the two time 
intervals, Tq, and Tngp. Tq is the a relaxation time defined 
by Fs{km,Ta) = e~^, where Fs{k,t) is the self-part of 
the density time correlation function for particle species 
1, and km = 27r is the first peak wavenumber of the 
static structure factor. Tngp is the time at which the non- 
Gaussian parameter [24] of the van Hove self-correlation 
function is maximized. In Fig. 2, we show Tq, and Tngp 
as function of the inverse temperature l/T. Ta — Tngp 
at T = 0.306, and Tq grows exponentially larger than 
Tngp with decreasing temperature at T < 0.306. This 
trend agrees with other simulation results of Lennard- 
Jones (LJ) systems [25] [26]. Next, we visualize the het- 
erogeneous dynamics of t^ and Tngp in the same manner 
presented in Ref.[6]. We calculate the displacement of 
each particle of species 1 in the time interval [to, to + 1], 
^rj{to,t) = rj{to + t)-rj{to) {j = 1, 2, iVi). In Fig.3, 
particles are drawn as spheres with radii 



a^(to,i) 



|Ar,(to,t)|' 
[Ar,(to,t)]2) 



(4) 



located at Rj{to,t) = \ [i'.j{to) + i'j{to + 1)] in both time 
intervals: [to, to + Ta] {t = t^) in 3(a), [^o, to + Tngp] {t ~ 
Tngp) in 3(b). The temperature is 0.253. a'j{to,t) > 1 
{a'j{to,t) < 1) means that the particle j moves more (less) 
than the mean value of the single-particle displacement. 
In Fig.3, the red (blue) spheres represent a'j{to,t) > 1 
{aj{to, t) < 1). We can see the large-scale heterogeneities 
in both Ta and Tngp. 

We can quantify the lifetimes of the heterogeneous dy- 
namics in both time intervals. We consider the local fiuc- 
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FIG. 3. Visualization of the heterogeneous dynamics of particle species 1. The temperature is 0.253. The time interval is 
[to,to + Ta] in (a), [to , to + Tngp] in (b). The radii of the spheres are \Arj{to,t)f/{[Arj{to,t)]^), and the centers are at |[^i(to) + 
rj{to + t)]. The red spheres and blue spheres represent \Arj{to,t)\ / {[Arj{to,t)]'^) > 1 and \Arj{to,t)\ / {[Arj{to,t)]^) < 1, 
respectively. 



tuations in the particle dynamics defined by 



SV{q,ta,t) 



Ni 

i=i 



ia^Jto,t) - 1) exp[-iq- • i?, (to, t)], (5) 



which is equal to the fluctuations in the diffusivity density 
defined in ref [6] and represents the local fluctuations in 
the particle dynamics in the time interval [to, to + t\. We 
use ST>{q,to,t) as SQkiq,to,t) in Eq.(2), and the time 
correlation function defined by 

Sviq, ts,t)^ {5V{q, t, + t, t)5V{-q, 0, i)), (6) 

corresponds to Fi^k{q_,ts,i)- Sv{(l,ts,i) represents the 
correlation of the particle dynamics between two time 



intervals [0,t] and [^s +t,ts + 2t] (see Fig.l). Thus, we 
can estimate the lifetime of the heterogeneous dynamics 
by examining the time decay of Sx>{q,ts,t). As the time 
separation increases, ST>iq,ts,t) with fixed t decays in 
the stretched exponential form, 



Sv{q,ts,t) 
Sv{q,0,t) 



exp 



Th{q,t) 



(7) 



where Th{q,t) is the wavenumber-dependent heteroge- 
neous dynamics lifetime, which corresponds to Ti^k{q,i) 
in Eq.(3). Figure 4 shows the time decay of ST>{q,ts,t) 
at t = Ta, q = 0.38 for various temperatures, q = 0.38 is 
the smallest wavenumber in our simulation. In Fig. 5, we 
show the wavenumber dependence of Th{q,t) aX t = Tq.. 
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FIG. 4. The time decay of S-D{q,ts,t) a.t t = t^, q = 0.38 for 
T — 0.772 - 0.253. g — 0.38 is the smallest wavenumber in 
our simulation. Temperature decreases going right. 



FIG. 5. The wavenumber dependence of Th{q,t) a.t t — 
for T = 0.772 — 0.253. Temperature decreases going up. The 
dotted line is Th{q,Ta) ~ q~^ 
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FIG. 6. The lifetime rhetero(i) for t = Tc, Tngp versus Ta- 
The Hne Thetero = 2.65r^ '^* is fitted for t = Ta, while the line 
Thetero = 3.75r°'^^ is fitted for t = Tngp. 

ThiQ, Ta) depends on q more weakly than the g-dependent 
relaxation time of the two-point density correlation func- 
tions and dramatically increases with decreasing temper- 
ature in a wide region of q (q ^ 0.38 — 19). Furthermore, 
we can see that Th{q,Ta) approaches Th{q,Ta) ~ q~^ 
at small wavenumbcrs. This suggests that the hetero- 
geneous dynamics migrate in space with a diffusion-like 
mechanism. These results of Th{q,Ta) are qualitatively 
the same as those of Th{q,Tngp)- 

We determined the lifetime of the heterogeneous dy- 



namics rhctcro(i) as Th{q, t) sd, q = 0.38, which is the time 
separation tg at which Sx){q,ts,t)/ Sxi{q,0,t) at q = 0.38 
equals e^^ in Fig.4. rhetero(i) increases dramatically 
with decreasing temperature. We plot rhotoro(i) versus 
Ta in Fig.6, which shows that 

'^hctcro (Ta) t1-08 and 
Thetero (Tngp) ^ T^'^^ . The difference between Thetero(Ta) 
and rhetero(7ngp) increases with decreasing temperature. 

At T = 0.253, rhetero(Ta) ^ 7.8Tq, and rhetero(Tngp) ^ 

1.4tq,. Therefore, Thetero (Ta) is considerably larger than 
Ta, while Thetero (Tngp) is Comparable to t^. The existence 
of a slower timescale in the heterogeneous dynamics is 
consistent with Ref.[20]. 

Finally, we examine the finite-size effect. To this end, 
we performed MD simulations using a larger system with 
Ni = N2 = 50,000 and L = 50ai, and compared our 
results with those of a larger system. No finite-size effect 
was detected in quantities such as Tq,, Tngp, or Thetero- 

In summary, we have investigated the heterogeneous 
dynamics in two different time intervals, Tq, and Tngp. We 
quantified the lifetimes of the heterogeneous dynamics in 
these two intervals, Thetero (to) and Thetero (Tngp), by cal- 
culating the time correlation function of the particle dy- 
namics. We found that the difference between Thetero (to) 
and Thetero(Tngp) iucrcascs with decreasing temperature. 
At low temperatures, Thetero (tq) is considerably larger 
than Ta, while Thetero (Tngp) remains comparable to Ta- 
Thus, we can conclude that the lifetime of the heteroge- 
neous dynamics depends strongly on the time interval. 
We also have examined the finite-size effect. No finite- 
size effect was detected in our study. 
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